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Abstract 

Background: The polymorphism rs6313 (T102C) has been associated with suicidal behavior in case-control and 
meta-analysis studies, but results and conclusions remain controversial. The aim of the present study was to 
examine the association between T102C with suicidal behavior in a case-control study and, to assess the combined 
evidence - this case-control study and available data from other related studies - we carried out a meta-analysis. 

Methods: We conducted a case-control study that included 161 patients with suicide attempts and 244 controls; 
we then performed a meta-analysis. The following models were evaluated in the meta-analysis: A) C allele vs T 
allele; B) T allele vs C allele; C) Caucasian population, D) Asian population, and E) suicide attempters with 
schizophrenia. 

Results: We found an association between attempted suicide and control participants for genotype 

(X2=6.28, p=0.04, df=2) and allele 0(2=6.17, p=0.01, df=1, OR 1.48 95% IC: 1.08-2.03) frequencies in the case-control 

study. The meta-analysis, comprising 23 association studies (including the present one), showed that the rs6313 

polymorphism is not associated with suicidal behavior for the following comparisonsT allele vs C allele (OR: 1.03; 

95% CI 0.93-1.13; p(Z)=0.44); C allele vs T allele: (OR:0.99; 95% CI: 0.90-1.08; p(Z)=0.22); Caucasians 

(OR:1.09; 95% CI: 0.96-1.23), and Asians (OR:0.96; 95% CI: 0.84-1.09). 

Conclusion: Our results showed association between the rs6313 (T102C) polymorphism and suicidal behavior in 
the case-control study. However, the meta-analysis showed no evidence of association. Therefore, more studies are 
necessary to determine conclusively an association between T102C and suicidal behavior. 
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Background 

Suicidal behavior is a major health problem worldwide. 
Several recent studies have been carried out that support a 
possible relationship between genetic factors and suicidal 
behavior [1,2]. Historically, evidence for the involvement 
of serotonin (5-HT) in suicide originated from findings of 
low 5-hydroxyindoleacetic acid concentration (5-HIIA) in 
cerebrospinal fluid (CSF) of depressed suicide attempters 
and in brain stems of completed suicides [3-5]. These 
studies provided evidence for altered serotonergic neural 
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transmission in the pathogenesis of suicidal behavior [6,7]. 
In consequence, genes pertaining to the serotonergic 
system have been proposed as candidates to establish bio- 
logical correlates between suicidal behavior and the sero- 
tonergic system [8-12]. 

Fourteen known serotonin receptor subtypes are 
involved in serotonin action [13]. One candidate gene in 
the study of suicidal behavior is the gene encoding the 
serotonin 2A receptor. There is evidence that the density 
of serotonin 2A receptors is upregulated in parietal cor- 
tical regions of depressed suicide victims and this in- 
crease has been suggested to be at least partly regulated 
genetically [14]. The 5HTR2A gene is located on 
chromosome 13 (13ql4-q21); it spans 20 kb and 
contains 3 exons; more than 200 SNPs have been identi- 
fied along the gene [13,15]. Mainly, the polymorphisms 
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A-1438G (rs6311), T102C (rs6313) and His452Tyr 
(rs6314) have been studied in relation to suicidal behav- 
ior. Moreover, a recent report has studied all SNPs in 
5HTR2A [16]. One of the most remarkable single nu- 
cleotide polymorphisms of the 5HTR2A gene is T102C 
(rs6313). Several genetic studies have shown an associ- 
ation between 5HTR2A and suicidal behavior [17,18]. 
Similarly, this common variant has been studied not 
only in relation to suicidal behavior, but also in associ- 
ation with other diseases such as schizophrenia [19-21], 
alcoholism [22], and depression [23]. 

To our knowledge, up to the present moment, there are 
no association studies of 5HTR2A and suicidal behavior in 
the Mexican population. Therefore, we examined the asso- 
ciation between 5HTR2A and suicidal behavior in this 
population. This work focuses on investigating the genetic 
basis of suicide attempters by assessing the rs6313 geno- 
type. Finally, we used the combined evidence to carry out 
a meta-analysis of all published data. 

Methods 

Case-control study 

A total of 161 patients were consecutively recruited 
from the outpatient service of the General Hospital of 
Comalcalco in the state of Tabasco, Mexico. These 
patients had attempted suicide between January and 
December 2011. In addition, 247 unrelated controls 
participated in this study. They were recruited from 
the Blood Donor Center of the General Hospital of 
Comalcalco and from the general population of the 
Comalcalco city area in the state of Tabasco, Mexico. 
Subjects were physically healthy on medical evaluation; 
none manifested psychiatric problems, as assessed in 
brief interviews by psychiatrists. To reduce ethnic vari- 
ation and stratification effects, only Mexican subjects 
descending from Mexican parents and grandparents 
participated in this study. However, no structured meth- 
ods were used to assess genetic heterogeneity. 

Ethics statement 

After they were given a verbal and written explanation of 
the research objectives, all subjects signed an informed 
consent to participate in the study; they did not receive 
any economical remuneration. This study complied with 
the principles convened in the Helsinki Declaration. In 
addition, this study was approved by the DAMC-UJAT 
Ethics and Research Committee (P.O.A. 20110237). 

Clinical evaluation 

Suicide attempt patients were evaluated by a trained 
psychiatrist or clinical psychologists with at least a mas- 
ter's degree level, using the Structured Clinical Interview 
for DSM-IV Axis-I and II diagnoses, in the Spanish ver- 
sion. We defined suicide attempt as a self-harm behavior 



with at least some intent to end one's life [24]. Subjects 
were excluded when self-injury behaviors were deter- 
mined to have no suicidal intention or ideation, based 
on the Assessment of suicidal intention: the Scale for 
suicide ideation in the Spanish version [25]. 

Genotype assays 

Genomic DNA was extracted from peripheral blood 
leukocytes using a modified version of the protocol by 
Lahiri [26]. The 5HTR2A T102C (rs6313) genotypes 
were analyzed in all patients using a polymerase chain 
reaction (PCR) end-point method. The final volume of 
the PCR reaction was 5 uL and consisted of 20 ng 
genomic DNA, 2.5 Fluorescence Labeling (FL) TaqMan 
Master Mix, and 2.5 FL 20x Assay. The amplification 
was performed in 96-well plates using the TaqMan Uni- 
versal Thermal Cycling Protocol. After the PCR end- 
point was reached, fluorescence intensity was measured 
with the 7500 Real-Time PCR system using SDS v2.1 
software (Applied Biosystems). All geno typing was per- 
formed blind to patient outcome. As a quality control in 
our genotyping analyses we used random blind duplicates. 

Statistical analysis 

Hardy- Weinberg equilibrium was tested using Pearson's 
goodness-of-fit chi-squared test. Chi-squared test or 
Fisher's Exact test was used to compare genotype and al- 
lele frequencies between groups. The power to detect 
associations given the sample size was analyzed by using 
the Quanto 1.2 software. The power of the analysis was 
0.43 (minor allele frequency= 0.3, type of model: domin- 
ant, effect size: 1.5). The level of significance was set at 
p=0.05. 

Meta-analysis 

This meta-analysis followed the Preferred Reporting Items 
for Systematic Reviews and Meta-Analyses (PRISMA) cri- 
teria [27]. 

Identification and selection of publications 

A literature search was performed that covered from 1997 
to 2012. Relevant publications were identified using the 
following search terms in Medline, Ebsco and Web of Sci- 
ence databases: "HTR2A and suicidal behavior", "HTR2A 
and suicide", "rs6313 and suicide", "5-hydroxytryptamine- 
receptor 2A" and "serotonin receptor 2A". These words 
were combined to retrieve the summaries. The search also 
implicated the review of the bibliography cited at the end 
of various research articles to identify additional papers 
not covered by the electronic search of abstracts. The pub- 
lications had to fulfill the following selection criteria: 
(1) to be studies on the relationship between rs6313 and 
suicidal behavior (attempted, ideation and completed sui- 
cide), (2) to be published in peer-reviewed journals, (3) to 



Gonzalez-Castro et al. BMC Psychiatry 2013, 13:25 
http://www.biomedcentral.com/1471-244X/13/25 



Page 3 of 10 



contain independent data, (4) to be case-control associ- 
ation studies in which the frequencies of three genotypes 
were clearly stated or could be calculated, (5) the use of 
healthy individuals as controls, and (6) to provide suffi- 
cient data to enable the estimation of an odds ratio with 
95% confidence interval. 

Two investigators (Gonzalez-Castro TB and Tovilla- 
Zarate CA) working independently screened each of the 
titles, abstracts and full texts to determine their inclu- 
sion. The results were compared and disagreements 
were resolved by consensus. 

Data extraction 

The same authors mentioned in the last paragraph care- 
fully extracted the information from all the included 
publications; they worked in an independent manner 
and in accordance with the inclusion criteria listed 
above. The following data were obtained from each of 
the studies: authors, year of publication, region, number 
of cases and controls, age, and ethnical origin of the par- 
ticipants. The outcomes of the meta-analysis were built 
by taking into consideration the following categories: a) 
exposed sick, b) exposed not-sick, c) not-exposed sick 
and d) not-exposed not-sick. The "sick" term refers to 
subjects exhibiting suicidal behavior and the "exposed" 
term to the allele of risk. 

Data analysis 

For the meta-analysis procedures, we used the EPIDAT 
3.1 program (http://dxsp.sergas.es). This software is freely 
available for epidemiologic analysis of tabulated data. Data 
were analyzed with the random-effects model following 
the reports in the literature [2,28,29]. Sample heterogen- 
eity was analyzed with the Dersimonian and Laird's Q test. 
Q test results were complemented with graphs to help the 
visualization of those studies favoring heterogeneity. The 
results of the meta-analysis are expressed as an odds ratio 
(OR). To address the problem of publication bias, funnel 
plots were calculated with the EPIDAT 3.1 software. This 
plotting standardizes the effect of each of the published 
studies on the vertical axis and its correspondent precision 
on the horizontal axis. Finally, a chi-squared (x 2 ) analysis 
was used to calculate the Hardy- Weinberg equilibrium to 
assess genotype distribution. 

Results 

Case-control study 

In the present study a total of 161 patients (61 males, 
100 females) were recruited; mean age was 25.5 (9.56) 
years old (range: 14-56 years of age). Control subjects 
consisted of 244 volunteers (82 males, 162 females); 
mean age was 33.1 (13.0) years old (range: 14-61 years 
of age). Both groups (cases p=0.06 and controls p=0.06) 
showed Hardy- Weinberg equilibrium for the analyzed 



genetic variability. From the 161 suicide attempters, 11 
(7%) exhibited the T102T genotype, 80 (50%) the T102C, 
and 70 (43%) the C102C type. On the other hand, in the 
control group, 9 individuals (4%) presented the T102T 
genotype, 100 (41%) the T102C type, and 135 (55%) the 
C102C type. We observed significant differences in 
genotype 0(2=6.28, p=0.04, df=2) and allele (x2=6.17, 
p=0.01, df=l, OR 1.48 95% IC: 1.08-2.03) frequencies be- 
tween patients and the control group. 

Meta-analysis study 

The selection process in this study is detailed in Figure 1. 
With regard to the literature search, a total of 397 studies 
were identified, but only 23 were used in this meta-analysis, 
including our case-control study [14,20-23,30-40]; this 
involved 2566 cases and 3989 controls (Table 1). In 
addition, 612 cases from a more recent meta-analysis and 
1129 controls were incorporated as well [41]. Overall, the 
meta-analysis included 13 studies in Caucasian populations, 
6 studies comprising Asian populations and 4 from other 
ethnic origins. 

We measured the Hardy- Weinberg equilibrium in all 
genotyped populations. Both patients and controls were 
in equilibrium (p>0.05), with the exception of the con- 
trol group in Du et al. (p=0.004), Turecki et al. (0.009) 
and Crawford et al. (0.01) [14,30,35], as well as the 
patients in Zhang et al. (p=0.04) [38]. When we explored 
all populations in a combined way we still encountered 



397 potentially relevant articles 
(Pudmed =1 57 Ebsco 1 1 4 Web of Science= 1 26) 



Excluded based on title 
n= 203 



49 references screened by abstracts 



Excluded 
n= 4 review articles 
n= 18 no outcomes of interest 



27 relevant studies retrieved for detailed 
evaluation 



Excluded n= 4 



23 studies in the meta-analysis 



Figure 1 Flow-chart showing the different phases of the 
meta-analysis. 



Table 1 Descriptive characteristics from 23 studies on the role of theT102C polymorphism of the 5HTR2A gene in suicidal behavior 



Study (Year) 


Sample 
Size n 
(Cases- 
Controls) 


Location 


Diagnosis 


Number 
of C 
alleles 


Gender (males/females) 
Cases Controls 


Mean age (years) 
Cases Controls 


Tan (2002) 


47/94 


Singapore 


SA with schizophrenia 


32 


47 


94 






Pooley (2003) 


1 29/329 


United Kingdom Caucasians 


SA 


154 


52/77 


138/191 


38 


38 


Chen (2001) 


471/523 


Hong Kong 


SA with schizophrenia 


193 


314 


308/215 






Saiz (2008) 


1 93/420 


Asturias in Northern Spain 


SA 


196 


70/1 23 


216/204 


35 


40 


Du (2000) 


42/131 


Ottowa 


SA with major depression 


99 


48/30 


61/70 


40 


36 


Wrsozek (2011) 


38/118 


Polish Caucasians 


SA with alcohol-dependent 


33 


26/12 


81/29 


>18 


>18 


Khait (2005) 


52/63 


European Caucasians 


SA with depression 


54 


30/24 


12/51 


37 


40 


Correa (2007) 


42/85 


Brazilian 


SA with schizophrenia 


83 


23/19 


45/40 






Vaquero-Lorenzo (2008) 


441/410 


Spanish Caucasians 


SA 


472 


161/280 


237-173 


38 


36 


Du (1999) 


24/31 


Budapest, Hungary 


SC with depression 


32 


24 


31 


45 


50 


Bondy (2000) 


131/125 


Caucasians Germany 


SA with major depression 


142 


91/40 


60/75 


'16 


47 


Zhang (2008) 


297/303 


Shangai 


SA with schizophrenia, anxiety and organic mental disorder 


229 


171/126 


217/112 


46 


43 


Preuss (2000) 


45/1 1 7 


Caucasians Germany 


SA with alcohol-dependence 


52 


45 


56/61 




47 


Preuss (2000) 


62/1 1 7 


Caucasians Germany 


SI with alcohol-dependence 


66 


62 


56/61 




47 


Ertugrul (2004) 


71/26 


Caucacians and African-americans 


SA with schizophrenia 


49 


26/12 


81/29 






Crawford (2000) 


78/131 


Caucasians Australian 


SC 


73 


68 


95 






Crawford (2000) 


68/95 


Caucasians Australian 


SI 


99 


78 


131 






Turecki (1999) 


56/1 26 


Montreal 


SC 


68 


56 


126 


36 


33 


Ono (2001) 


151/163 


Japanese 


SC 


157 


1 06/45 


108/55 


47 


47 


Zhang (1997) 


15/150 


Japanese 


SA with mood disorders 


21 


15 


65/85 


33 


37 


Tsai (1999) 


61/96 


Taipei 


SA with mood disorders 


49 


61 


96 




71 


Arias (2001) 


33/164 


Spanish 


SA 


42 


33 


83/71 


55 


41 


Tovilla (2012) 


161/244 


Mexican 


SA 


220 


61/100 


82/162 


25 


33 



SA: Suicide attempters; SC: Suicide completers; SI: Suicide ideation. 
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Study (Year)n- 
Zhang (1997) 330- 
Du (1999) 110 - 

Tsai (1999) 314- 

Turecki (1999) 250 - 
Bondy (2000) 512 - 

Du (2000) 418 - 
Preuss (2000) 324 - 
Preuss (2000) 358 - 
Crawford (2000) 326 - 
Crawford (2000) 418- 
Arias (2001) 394- 
Chen (2001) 1674- 

Tan (2002) 282 - 

Pooley(2003) 142 - 
Ertugrul (2003) 142 - 
Vadim (2005) 234 - 
Correa (2007) 258 - 
Zhang (2007) 1084- 
Vaquero- Lorenzo (2008) 1702- 
Wrzosek(2D11) 277 - 
Tovilla (2012) 810- 



OR Cl(95.0%) 
2.63(1.16,5.93) 
0.88(0.39,1.98) 
1.16(0.73,1.86) 
5.40(3.11,9.34) 
0.96(0.67,1.36) 
1.81(1.21,2.72) 
1.15(0.70,1.88) 
0.95(0.61,1.48) 
0.64(0.41,1.01) 
1.81(1.21,2.72) 
1.31(0.76,2.28) 
0.87(0.70,1.08) 
1.04(0.62,1.77) 
0.90(0.67,1.21) 
1.25(0.63,2.47) 
1.21(0.72,2.02) 
0.86(0.51,1.46) 
0.91(0.71,1.15) 
1.22(1.00,1.47) 
0.89(0.49,1.62) 
0.68(0.50,0.94) 



Bondy(2000) 512- 

Du (2000) 418 - 
Preuss (2000) 324 — 
Preuss (2000) 358— 
Crawford (2000) 326 - 
Crawford (2000) 418 - 
Arias (2001) 304 - 
Chen (2001) 1674- 

Tan (2002) 282 - 

Pooley (2003) 916 - 
Ertugrul(2003) 142- 
Vadim (2005) 234 - 
Correa (2007) 258 - 
Zhang (2007) 1084- 
Vaquero- Lorenzo (2008) 1702 — 
Wrzosek (2011)277- 
Tovilla (2012) 810 - 



6(0.67,1. 

1(1.21,2.72) 

5(0.70,1. 

5(0.61,1.48) 

4(0.41,1.01) 

1(1.21,2.72) 

1(0.76,2.2; 

7(0.70,1.01 

4(0.62,1.77) 

0(0.67,1.21) 

5(0.63,2.47) 

1(0.72,2.02) 

6(0.51,1.46) 

1(0.71,1.15) 

2(1.00,1.47) 

9(0.49,1. 



Random effects 1 1133 - 



Random effects11133 - 



2 3 4 5 6789 



I I I I I I I 



Study (year) n — 
Turecki (1999) 250 - 
Bondy (2000) - 

Du (2000) 418- 

Preuss (2000) 324 - 
Preuss (2000) 358— 
Crawford (2000) 326— 
Crawford (2000) 418— 
Arias (2001) 394- 
Ertugrul (2003) 142- 
Pooley (2003) 916 - 
Vadim (2005)n 324 — 
Vaquero- Lorenzo (2008) 1702 — 
Saiz(2008) 1226 — 
Wrzosek (2011)277 — 



Random effects 7497- 




Weight (%) OR Cl(95%) 
6.15 



Zhang (1997) 330 — 



□u (1999) 110- 



Tsai (1999) 314 



Chen (2001) 1674- 



Zhang (2007) 1 084 - 



Random effects 4422 - 



0.8 1 1.2 1.4 1.6 1.8 2 2.2 2.4 



B 




Welght{%) OR CI (95%) 



4.21 2.63(1.16,5.9! 



4.26 0.88(0.39,1.98) 



1083 1.16(0.73,1.86) 



27.83 0.87(0.70,1.08) 



1877 1.16(0.85,1 



9.01 1.04(0.62,1.77) 



25.05 0.91(0.71,1.15) 



C D 

Figure 2 Odds ratios and forest plots of all models in overall studies. A) C allele vs T allele; B) T allele vs C allele; C) Caucasian population, 
and D) Asian population. 
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them in equilibrium with p=0.19 and p=0.92 for controls 
and cases, respectively. 

Figure 2 shows the pooled OR derived from all stud- 
ies indicating a non-significant association of the C 
allele in the T102C polymorphism with suicidal behav- 
ior (Random effects: OR:1.14; 95% CI 0.96-1.35; p(Z) 
=0.07). We observed heterogeneity in all studies 
(Q=75.03 df=20; p=<0.0001). The Egger's test provided 



no evidence of publication bias (t=1.37, df=19; p=0.18) 
(Figure 3A). Therefore, we carried out a second ana- 
lysis, which excluded studies based on the heterogen- 
eity curve and sensitivity analysis. As a result, the 
reports by (Crawford et al. [35], Zhang et al. [38], Du 
et al. [30], Turecki et al. [14] and ours were excluded. 
However, we could not find an association either (OR: 
1.03; 95% CI 0.93-1.13; p(Z)=0.44). 
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Table 2 Meta-analysis of case-control studies on the role of the T and C alleles in the T102C polymorphism of the 
5HTR2A gene in individuals with suicidal behavior 



Rpfpronr 


Number of T alleles 


Number of C alleles 


Odds Ratio (95%IV) 






V_OI 1 IT \Jl3 




V-UI 1 LI UI 5 


T Allele 


C Allele 


CoTr^a 


91 


A 1 


8 -I 


43 


1 1 4 (0 ftR - 1 Q31 


1 14 (C\fc,R - 1 Q31 


V V IZUifcrrx 


24 


8 / 


33 


1 33 


1 . 1 1 \\J.U 1 Z..UU) 


1 1 1 ( n fti 1 001 

1 . 1 1 1 U.D 1 Z.UUJ 




287 


303 


229 


265 


1 DQ fn RfS - 1 3Q1 

1 .U:7 ^U.OU 1 .Dzj) 


1 OQ ( DRfi - 1 3Q1 


vao, i m 


54 


69 


54 


57 


0 9,7 lT)4Q - 1 381 

U.Oi ^U/TI7 1 JO/ 


OR? f D4Q - 1 3R1 

U.OZ ^ U/TI7 1 JO/ 


Du 


5 7 


1 34 


99 


1 28 


0 ^ (0 3ft - 0 89^41 
U.JJ 1,U.jU U^ZDH-^ 


1 1 3 (0 ^0 - 9 ^31 

1 . 1 J \.U.JU Z..DD) 


VdL|Ucl U LUl ci 1ZU 


4 1 0 


422 


4 72 


398 


0 R1 fO fl- 0 QQ1 
U.o 1 [U.U/ \J.yy) 


0 R1 I 1 0 rS7- 0 QQ1 
U.o l ^ u.u/ u.yyj 


Du 


1 6 


1 9 


32 


43 


1 1 3f 0 ^0 - 1 ^31 
1 . 1 j\ U.JU Z.J JJ 


n^w n R9i 

U.JJ I. U.JO U.OZJ 


Rj~inH\/ 
□ui luy 


120 


1 12 


142 


138 


~] f)A ( n 73- 1 4.71 


1 04 ffl 73 - 1 471 

1 .vrr ^U./ J 1 / / 


Zhang 


Q 


1 59 


21 


1 4 1 


1 HQ fflRfi - 1 3Q1 
1 .U:? I^U.OO 1 J7j 


1 OQ (f) Rft - 1 3Q1 




24 


1 41 


'■12 


1 87 


0 7q / ri4^ - 1 301 

U.I J I. U.tO 1 .DU) 


0 3R ( n ifi - n Rm 

U.JO ^ U. 1 D U.OJ) 


Tsai 


/3 


1 22 


49 


70 


n or / pi r-2 _ 1 
U.OJ ^ U.JJ 1 .DU) 


0 R^f 0^3-1 3ft1 
U.OJl, U.JJ 1 .joy 


Turccki 


40 


108 


68 


34 


ms (o m - n 391 

U. 1 O \U. 1 U U.JZJ 


n is c n m - n 391 

u. 10 ^ u. 1 u u.jzy 




38 


1 0 / 


52 


1 27 


n r/s f n s ^ - 1 411 

U.OU ^ U.JJ LtIJ 


n c n 3A - nj39i 

U JJ V U.JU U.OZ/ 


Preuss^ 


58 


1 0 / 


66 


1 27 


1 04 ! 0 67 - If,]) 
1 .Lr+ ^ U.O/ 1 .O 1 J 


1 04 fO ft7 1 ft1 ) 
I .Ut \U.O/ I .U I J 


d VVI (J[ U 


63 


68 


73 


122 


1 .54 ( 0.98 - 2.42) 


1 .54 ( 0.98 - 2.42) 


Crawford 


57 


134 


99 


128 


0.55 ( 0.36 - 0.82) 


0.55 ( 0.36 - 0.82) 


Ertugrul 


39 


27 


49 


27 


0.79 (0.40- 1,56) 


0.79 (0.40 - 1.56) 


Tan 


62 


126 


32 


62 


0.95 (0.56- 1.60) 


0.95 (0.56 - 1 .60) 


Pooley 


104 


249 


154 


409 


1.10 ( 0.82 - 1.48) 


1.10 (0.82 - 1.48) 


Chen 


435 


694 


193 


352 


1.14 (0.92- 1.41) 


1.14 (0.92 - 1.41) 


Tovilla 


102 


118 


220 


370 


1.45 (1.06- 1.98) 


1 .45 (1 .06 - 1 .98) 


Ono 


145 


169 


157 


157 






Saiz 


190 


378 


196 


462 






Random effects 


2498 


3894 


2615 


3980 


0.87 ( 0.73 - 1.03) 


0.87(0.73 - 1 .03) 



Given that our results in the association study sug- 
gested an association with the T allele, we performed a 
meta-analysis using the allelic model for this allele 
(Table 2). We did not encounter any association between 
the T allele and suicidal behavior (Random effects; 
OR: 0.87; 95% CI 0.73-1.03; p(Z)=0.07). We also found 
heterogeneity in all studies (Q=75.03 df=20; p=<0.0001). 
The Egger's test indicated no evidence of publication 
bias (t=-1.37, df=19; p=0.18; Figure 3B). A second ana- 
lysis was performed, which only included studies based on 
the heterogeneity curve and sensitivity analysis (Random 
effects: OR: 0.99; 95% CI: 0.90-1.08; p(Z)=0.22). The Egger's 
test gave no evidence of publication bias (t=-0.1726, df=12; 
p=0.86). 

In addition, we performed an analysis in two popula- 
tion groups: Caucasian (Fixed effects; OR: 1.14; 95% 
CI: 1.03-1.25, Random effects OR: 1.21; 95% CI: 0.97- 
1.57 with heterogeneity; Egger's test Figure 3C, OR: 1.09; 
95% CI: 0.96-1.23 without heterogeneity; Table 3) and 
Asian (OR: 0.96; 95% CI: 0.84-1.09 without heterogen- 
eity; Egger's test Figure 3D). Again in this case, we could 
not find any association with suicidal behavior. Finally, 



in an analysis including only studies of suicide attempt 
patients with schizophrenia, the results were also nega- 
tive (Random effects; OR: 0.91; 95% CI: 0.79-1.06 with- 
out heterogeneity). 

Discussion 

In this study, we explored the association of the T102C 
(rs6313) polymorphism of the 5-HTR2A gene with sui- 
cidal behavior. First, a case-control study was con- 
ducted. Subsequently, we performed a meta-analysis to 
assess the current evidence of association between 



Table 3 Results of Caucasian and Asian populations with 
suicidal behavior 





OR 


CI (95%) 


Overall studies with heterogeneity 


1.14 


(0.96,1 .35) 


Overall studies without heterogeneity 


1.03 


(0.93,1.13) 


Caucasians with heterogeneity 


28.3704 


(16.91,47.57) 


Caucasians without heterogeneity 


1.09 


(0.96,1 .23) 


Asians with heterogeneity 


1.02 


(0.86,1 .22) 


Asians without heterogeneity 


0.96 


(0.84,1.09) 
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rs6313 and suicidal behavior. In the first part, we found 
an association between rs6313 and suicidal behavioral in 
a Mexican sample. To our knowledge, this is the first 
study addressing the genetic association between the 
T102C allele and suicidal behavior in a Mexican 
population. 

The genetic analysis revealed a slight preponderance of 
the T allele in the cases group. In 2000, Du et al. described 
an association between the 102C allele and ideation 
suicide [30]. Since 2000, growing interest on this issue 
prompted a variety of studies in Caucasian and Asian 
populations dealing with 5HTR2A and suicidal behavior. 
Although many studies have reported this association 
[21,22,33,42], there are other reports that have not 
encountered an association of either the T or C allele in 
rs6313 of 5HTR2A [20,23,31,32,35,38,43] with suicidal be- 
havior. The slight preponderance of the T allele is in agree- 
ment with other reports in the literature. In 1999, Tsai 
et al. [34] reported that the allele frequency of 102T is 
higher than that of 102C in both patient and control 
groups. Recently, Saiz et al. [36] found an excess of the 
102T allele in a group of 193 patients and a control group 
of 420 individuals when comparing non-impulsive suicide 
attempts with impulsive suicide attempts. This evidence 
suggests that the rs6313 polymorphism of the 5-HTR2A 
gene may be involved in suicidal behavior. However, no 
conclusive outcomes have yet been attained. In agreement, 
our results in the Mexican population suggest that allele 
102T of 5HTR2A could be associated with suicidal 
behavior. 

There are several possible explanations for the discrep- 
ancies in the results regarding this association. First, 
there are differences in the populations of patients. For 
example, Du et al. [30] and other authors [13,23] have 
studied patients with major depression and suicidal be- 
havior; Wrzosek et al. studied patients with attempted 
suicides and alcohol dependence [22]; other studies have 
analyzed suicidal behavior in patients with schizophrenia 
[20], and finally, there are studies such as ours that have 
analyzed only suicide attempters [8,31,36]. Second, there 
are differences in genetic heterogeneity. Several authors 
have shown that allele frequencies of the T102C poly- 
morphism may be ethnic-dependent [34,36,44]. How- 
ever, our sample was made of a relatively homogenous 
Mexican population, because only patients of the state 
of Tabasco with parents and grandfathers born in 
Tabasco were included in this study. In consequence, we 
considered necessary to carry out a meta-analysis of all 
the published evidence concerning rs6313 and suicidal 
behavior. 

In the current meta-analysis we utilized the allelic model 
and ethnicity groups. There are two previous meta- 
analyses -the last including 25 studies- that have explored 
the association between rs6313 and suicidal behavior 



[41,45]. Both meta-analyses reported no association be- 
tween rs6313 and suicidal behavior. However, several 
additional case-control studies concerning this relation- 
ship and using larger samples have been recentiy reported, 
but they are not included in the above-mentioned 
meta-analyses. Therefore, to get a more comprehensive 
estimation we included our case-control study. 

We performed a meta-analysis for the C (OR: 1.14; 95% 
CI 0.96-1.35) and T (OR: 0.99; 95% CI: 0.90-1.08) alleles 
separately; however, no association was observed. The 
results suggest no association between rs6313 (T102C) 
and suicidal behavior. One possibility for this lack of asso- 
ciation may lie in the different criteria used in the selec- 
tion of patients. In addition, samples sizes of several of the 
studies included in our meta-analysis are in the low range 
compared with genetic studies for other diseases [7,46]. 

There is the possibility that the C allele may be regarded 
as a factor risk in the Caucasian population, because in a 
first analysis we observed an association when we 
considered a fixed effects with heterogeneity. However, we 
adopted the random effects, which accounts for additional 
sources of inter-study variation when heterogeneity exists. 
This model is more conservative than the fixed effects, 
since the latter assumes the same true genetic effects, 
whereas the former assumes normally distributed 
effects and parametrizes inter-studies variation [41]. 
When we performed a second analysis, in which the 
studies giving rise to heterogeneity were discarded, no 
association was observed between rs6313 and suicidal 
behavior in a Caucasian population. Therefore, future 
studies on suicidal behavior comprising larger samples 
are important to determine this association. Similarly, 
the accuracy of the phenotype definition in the associ- 
ation studies is a relevant issue [41]. However, it is ne- 
cessary to take into consideration other aspects such as 
different inheritance patterns and gene X environment 
interaction [16]. 

Finally, we recognize some limitations in the present 
study; 1) the sample size of the case-control study is small 
and may not have sufficient power to detect an association 
between suicidal behavior and small effects polymorphism. 
In relation to our meta-analysis, there are also some con- 
siderations to bear in mind: first; included studies were 
limited to published articles; second, we did not perform 
sub-analyses including attempted or completed suicide; 
third, we did not carry out a subgroup analysis based on 
gender, and fourth, we did not report psychopathologies 
related to suicide attempters. 

Conclusions 

In conclusion, we encountered a significant association 
between the T102C polymorphism and suicidal behavior 
in a Mexican population. However, this result was not 
observed in the meta-analysis. Nonetheless, our meta- 
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analysis provides pooled data on a substantial number of 
cases and controls that may provide a better understand- 
ing in the association between rs6313 and suicidal be- 
havior. However, more comprehensive studies and larger 
samples are necessary to determine conclusively the 
presence of this association. 
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